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In gauge systems coupled to matter, the static potential flattens out at a scale where the confining string breaks 
by formation of a dynamical pair of particles. Surprisingly, such a breaking is invisible in Wilson loops even when 
the inter-charge separation is much larger than the flattening scale. Observing string breaking also requires using 
different operators. A known string mechanism provides us with a simple explanation, leading to the area law 
for large Wilson loops, as observed in most gauge models coupled to whatever kind of matter. It is also pointed 
out that in a simple 3D Z-2 gauge-Higgs model, once reformulated in terms of Fortuin Kasteleyn clusters, this 
peculiar behaviour of the Wilson operator can be ascribed to non-trivial linking of a percolating cluster to the 
Wilson loop. Some numerical tests on this model are also presented. 



1. INTRODUCTION 

In gauge theories coupled to matter (quarks or 
Higgs fields in the fundamental representation) 
the string breaking is the phenomenon of static 
potential flattening at large distances due to the 
screening of the sources produced by pair cre- 
ation. A similar screening is expected also in pure 
Yang-Mills theory for sources in the adjoint rep- 
resentation, where the role of charged matter is 
played by the gluons (adjoint string). 

The coupled systems of this kind studied up to 
now reveal a surprising, general phenomenon: the 
string breaking is invisible in the Wilson loops, 
namely, the area law continues to prevail in full 
QCD [§,[§] and in pure Yang Mills theory with 
adjoint sources || even at distances where the 
static charges are completely screened. In other 
terms there is no possibility to extract the true 
large distance behaviour of the static potential us- 
ing only the Wilson loops. On the contrary, neat 
descriptions of potential flattening have been ob- 
tained in studies where the basis of the operators 
has been enlarged in order to gain a better ground 
state overlap. In this way it has been observed the 
breaking of the confining string between colour 
sources in fundamental representation in Higgs 
models Q, in QCD || and the breaking of the 
adjoint string ||. 



A possible explanation for this unexpected 
finding is that the available range of the Wil- 
son loops of present studies is not wide enough 
and that for larger sizes the overlap of the Wil- 
son operators to the broken string state could be- 
come visible. This however would imply that the 
presently observed area law obeyed by the Wil- 
son loops, which has been tested for a wide range 
even in QCD M, should be broken at a larger 
scale of doubtful physical significance. 

There is another explanation, based on the 
effective string picture of confinement and sup- 
ported by a well known string mechanism [j7|, 
which points in the opposite direction Ac- 
cording to such a picture, the area law describes 
the asymptotic, universal behaviour of the Wil- 
son loops in most (perhaps all) gauge theories in 
the confined phase, coupled to whatever kind of 
matter. Any observed non- vanishing overlap of 
the Wilson operator to the broken string state 
should be a decreasing function of the loop size, 
going eventually to zero in the IR limit. If this 
prediction will pass the numerical tests of the fu- 
ture lattice simulations, one can envisage using 
the asymptotic area law of the Wilson loops for 
an unambiguous definition of the confining phase 
and string tension even when the static potential 
flattens because of sea quarks or other dynamical 
matter. This will be discussed in section 2. 
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There is one case, at least, where the area law 
of the Wilson loops can be proved without ap- 
pealing to any effective theory: In the 3-D Zi 
gauge Higgs model it can be shown that the very 
confining mechanism which is operating in pure 
Zi gauge theory survives to the addition of Zi 
charged matter. In particular, for a simple topo- 
logical reason the Wilson loop continues to obey 
the area law. This will be discussed in section 3, 
where a very efficient numerical method is used 
to test the behaviour of large Wilson loops. 

2. EFFECTIVE STRING BREAKING 

The lack of visible effects of string breaking in 
the Wilson loop and its manifestation in other 
operators have a simple explanation in terms of 
the effective string picture of the confining phase 
of any gauge theory. According to such a de- 
scription, the expectation value of a rectangular 
Wilson loop can be represented as 
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where F Q is the free energy of a 2D model describ- 
ing the normal modes of vibration of the string 
world-sheet bounded by the rectangle RxT. This 
assumption yields, in the IR limit, 
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where q = exp(— 27rT/i?), a is the string tension 
and d the space-time dimension. Adding dynam- 
ical matter fields induces the formation of holes 
of any size in the world-sheet, hence eq. (|lj) is re- 
placed by a loop expansion: 
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where the loops inside the Wilson rectangle rep- 
resent pair created particles. What we can say 
about the sum of this series? A simple solvable 



matrix model, describing a free string with holes 
of any size, suggests that there are two phases: 

In the normal phase this sum is dominated by 
configurations with small holes. Their size does 
not depend on that of the Wilson loop: Larger 
Wilson loops have more holes. These holes do 
not influence the area law and their effect can be 
absorbed into a renormalization of the string ten- 
sion, according to a well known string mechanism 

§■ 

In the other (tearing) phase the mean size of 
the holes increases with that of the Wilson loop. 
In this case asymptotic screening should be visible 
also in the Wilson loops. 

The observed poor overlap of the Wilson oper- 
ator with the broken string state suggests that in 
all the known cases the confining string belongs 
to the normal phase. In this phase large Wilson 
loops behave exactly like in the pure Yang Mills 
theory. There are however finite size effects which 
are worth mentioning. As a matter of fact, differ- 
ent phases are completely separated only in the 
thermodynamic limit. The expectation value of 
finite Wilson loops in the normal phase receives a 
non- vanishing contribution of configurations typi- 
cal of the tearing phase, of course. Thus the over- 
lap of a finite Wilson loop to the broken string 
state cannot be exactly zero, but should be a de- 
creasing function of its size, going to zero in the 
infinite size limit. In order to gain an intuitive un- 
derstanding of what is going on, imagine to cut a 
rectangular Wilson loop R x T at a given "time" 
(the dashed line in the following drawing) 



o 


O CD O 


o o o 


o 


o 


o 


o cb o 


o o o 


o 


o 


O 


O i O 


O O O 


O 


o 


O CD O 


o o o o 


o 


o 


O 


oc^o 


O o 





o 



R 



T 



The above naive picture of the string world-sheet 
riddled with holes suggests that the Wilson oper- 
ator has an important overlap with a multi-meson 
state \n) where the mean number n of particles is 
a growing function of the distance R between the 
static sources. 
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The above reasoning does not imply that there 
is no flattening in the static potential. The 
observed screening of sources is a mixing phe- 
nomenon. It is simply due to the possibility, in 
presence of charged matter, to construct other 
operators with a sufficient overlap to the ground 
state even for large separations between the static 
source [f|, ||. Applying the loop expansion (||) 
to operators made of disconnected pieces leads 
to a contribution with a large overlap with the 
broken string state. This very argument also ex- 
plains why || in the QCD at finite temperature 
the string breaking is neatly visible already in the 
Polyakov correlator gj . 

3. 3D Z 2 GAUGE-HIGGS MODEL 

The action of a 3D Z 2 gauge theory coupled to 
a charged matter field can be written as 

S(J3 G , Pi) =PiJ2 +I 3 g/Z u o (4) 

(ij) plaq. 

where both the link variable f/y = Ui and the 
matter field Oi take values ±1 and Uq — Yiieu Ue- 
This model is self-dual: its partition function 

Z(0c,(3 I )= e ~ S(M) ( 5 ) 

{<7i=±l, U e =±l) 

fulfils, in the thermodynamic limit, the functional 
equation 

Z{Pg,Pi) = (sinh2/3 G sinh2/3 J )i^(/3 / J G ) (6) 

with $ = — \ log(tanh (3) . 

The phase diagram of this model has been stud- 
ied long ago Jlo| . There is an unconfined region 
surrounded by lines of phase transitions toward 
the Higgs phase and its dual. These lines are sec- 
ond order until they are near each other and the 
self-dual line, where first order transition occurs. 
We want to show that large Wilson loops in the 
dual Higgs phase obey the area law. 

Using the method of Fortuin and Kasteleyn 
(FK) we can map this system into a percola- 
tion model, rewriting Z as 

Z = e-P' N J2 e ^ G U ° E «" (G) 2 C(G) (7) 

Wil G 



with v = e 2 ' 1 — 1, G denotes a subgraph of the 
lattice made of n(G) links, called active bonds, 
and c(G) is the number of its connected compo- 
nents, called FK clusters. The only difference be- 
tween the pure gauge theory and the one coupled 
to matter is that in the former the G subgraphs 
are arbitrary, while in the latter the allowed G's 
are subjected to a constraint: the number of neg- 
ative links (Uij = — 1) along every circuit made 
with active bonds must be even, as it is easily 
checked. Put differently, no frustration is allowed 
in G, hence no Z 2 magnetic flux can pass through 
the circuits of G. This can be rephrased by saying 
that the FK clusters behave as pieces of super- 
conducting matter. On the other hand, insert- 
ing a Wilson loop W(R,T) in the vacuum cor- 
responds to creating an unit of Z 2 flux in the 
dual version. Thus, owing to the above super- 
conducting property, no FK cluster can be linked 
with it. Actually, using the methods of ref . Jl2| , 
one easily gets the following exact identity 

(W(R,T)) Pg!0i = {w(R,T))^ Jg (8) 

where w(R,T) is a projector on the G's defined 
as follows 

w(R T) — I ® ^ some FK cluster is linked to W 
I 1 if no FK cluster is linked to W. 

Thus the Wilson loop is a sort of cluster counter. 
It should be mentioned the surprising analogy of 
the circuts of the FK clusters with the center vor- 
tices |ll|,[[l4j]; note however that they have little 
to do with the thin vortices associated to the sign 
of the plaquettes; they are rather describing new 
degrees of freedom related to the gauge variables 
through a highly non-local duality transformation 
and live on links of the dual lattice. They play a 
relevant role in providing us with a very efficient 
disordering mechanism which leads to an area law 
for Wilson loops. 

For a large Wilson loop, finite clusters can link 
with it only along the loop perimeter. There- 
fore they contribute only to the perimeter term 
of eq.(|^). However, owing to the duality trans- 
formation involved in eq.(||), the Wilson loop in 
the dual Higgs phase is dominated by an infi- 
nite, percolating, FK cluster characterizing the 
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Higgs phase. It is precisely this percolating clus- 
ter which leads to the area law for large Wilson 
loops, as we wanted to show. 
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Figure 1. Universal shape effects in Wilson loops. 

Thus, in the dual of the Higgs phase, even if the 
ground state potential necesssarily flattens out, 
the area law behaviour of large Wilson loops is 
inavoidable. As a consequence we can give a pre- 
cise definition of confinement phase: it is the one 
where the area law holds true. In that phase the 
string tension a has an unambiguous definiton as 
the coefficient of the area term of eq.(^|). 

As a check of these results, we have estimated 
a in this gauge-Higgs model in a lattice of size 
L 3 = 40 3 at /3 G = 0.75245 and fa = 0.16683. We 
measured all the squared Wilson loops (W(R, R)) 
with 10 < R < 20 applying a very powerful 
algorithm already used in the pure gauge the- 
ory dD, |l6) and based on eq.((|). We found 
a = 0.00828(26). 

It also interesting to uncover the universal 
shape effects produced by the quantum fluctu- 
ations of th string. Actually the quantity 



R{nM= {W[L + n S L : n)) e-^ 



(9) 



(W(L,L)) 

according to the asymptotic form (^), is only a 
(known) function of the ratio does not con- 
tain any adjustable parameters, nor it depends 
on the gauge group or on the nature of the mat- 
ter fields. The function R(j;) is plotted in Fig.l 
(dotted line) along with the data of the present 
study and those of pure Z2 gauge theory taken 
form ref. Jl^ |. 
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